Topological Quantum Phase Transitions 
in Topological Superconductors 
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In this paper we show that BF topological superconductors (insulators) exibit phase transitions 
between different topologically ordered phases characterized by different ground state degeneracy 
on manifold with non-trivial topology. These phase transitions are induced by the condensation (or 
lack of) of topological defects. We concentrate on the (2+l)-dimensional case where the BF model 
reduce to a mixed Chern-Simons term and we show that the superconducting phase has a ground 
state degeneracy k and not . When the symmetry is 17(1) x U{1), namely when both gauge 
fields are compact, this model is not equivalent to the sum of two Chern-Simons term with opposite 
chirality, even if naively diagonalizable. This is due to the fact that U{1) symmetry requires an 
ultraviolet regularization that make the diagonalization impossible. This can be clearly seen using 
a lattice regularization, where the gauge fields become angular variables. Moreover we will show 
that the phase in which both gauge fields are compact is not allowed dynamically. 

PACS numbers: 74.20. Mn. 03.65.Fd, 05.30.Pr 



Quantum phase transitions describe changes in the 
entanglement pattern of the complex-valued quantum 
ground state wave function. The universality classes of 
these macroscopic quantum ground states define the cor- 
responding quantum orders When there is a gap in 
the spectrum, the quantum ordered ground state is called 
topologically ordered its remarkable hallmark is a 
ground state degeneracy depending only on the topology 
of the underlying space. 

The best known example of topological order is given 
by Laughlin's quantum incompressible fluids [3] describ- 
ing the ground states responsible for the quantum Hall 
effect [4j]. 

In [sl we proposed a superconductivity mechanism 
which is based on a topologically ordered ground state 
rather than on the usual Landau mechanism of spon- 
taneous symmetry breaking. Topologically ordered su- 
perconductors have a long-distance hydrodynamic action 
which can be entirely formulated in terms of generalized 
compact gauge fields, the dominant term being the topo- 
logical BF action. 

In this paper we will show that BF models exhibit 
phase transitions between different topologically ordered 
phases. These phase transitions are induced by the con- 
densation (or lack of) of topological defects. These de- 
fects are present due to the compactness of the original 
gauge model. Different topological order will be char- 
acterized by different ground state degeneracy. We will 
concentrate on the (2+l)-dimcnsional case where the BF 



model reduce to a mixed Chern-Simons term. We want to 
point out that when the symmetry is [/(l)xJ7(l), namely 
when the to gauge fields are compact, this model is not 
equivalent to the sum of two Chern-Simons term with 
opposite chirality, even if naively diagonalizable. This is 
due to the fact that U{\) symmetry requires an ultra- 
violet regularization. This can be typically treated by 
obtaining the U{\) group from spontaneous symmetry 
breaking of a non-Abelian group or by formulating the 
gauge theory on a lattice, in which case, the gauge fields 
are angular variables. In both regularization methods the 
ultraviolet cutoffs actually prevents the naive diagonal- 
ization of the mixed Chern-Simons term, in one case be- 
cause of the non-quadratic character of the action and in 
the other because of the angular character of both gauge 
fields. The diagonalization is, instead, permitted when 
only one of the gauge field is compact or when both are 
non-compact, i.e. for gauge groups U{1) x R, Rx U{1) 
OT R X R. As we have shown in ■()] and confirmed in 
the proper continuum formulation 7] these are actually 
the only cases that can occur due to the condensation 
of topological defects, the U{1) x U{1) case being dy- 
namically excluded. The three corresponding phases are 
superconducting, superinsulator and metallic phases. 

BF theories are topological theories that can be de- 
fined on manifolds M^+i of any dimension (here d is the 
number of spatial dimensions) and play a crucial role in 
models of two-dimensional gravity [8]. In [6*1 we have 
shown that the BF term also plays a crucial role in the 
physics of Josephson junction arrays. 

The BF term Q is the wedge product of a p-form B 
and the curvature dA of a (d-p) form A: 
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where fc is a dimensionless coupling constant. This ac- 
tion has a generahzed Abehan gauge symmetry under the 
transformation 



B 



where rj and ^ are a closed p and a closed (d-p) form, 
respectively. 

The degeneracy of the ground state of the BF theory 
on a manifold with non-trivial topology was proven in 
Consider the model ([T|) with fc = ^ on a manifold 
Md X i?i , with Md a compact, path-connected , orientable 
d-dimensional manifold without boundaries. The degen- 
eracy of the ground state is expressed in terms of the 
intersection matrix Mmn [l2 with m,n — l....Np and 
Np the rank of the matrix, between p-cycles and (d-p)- 
cycles. Np corresponds to the number of generators of 
the two homology groups Hp{Md) and Hd-p{Md) and is 
essentially the number of non-trivial cycles on the man- 
ifold Md- The degeneracy of the ground state is given 
by \kik2M\^'' , where M is the integer- valued determi- 
nant of the linking matrix. In our case p = (d-1) and the 
degeneracy reduces to 



fci X k2M 



(2) 



In this paper we consider the special case of (2+1) 
dimensions (d=2). In this case also B becomes a 1-form 
and, correspondingly the BF term reduces to a mixed 
Chern-Simons term with action 
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Ax A dBx . 
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In this case the degeneracy on a manifold of genus g is 
{k\ X ^2)^^, and thus (fci x k-if' on a torus 

In the application to superconductivity, the conserved 
current j\ = *dBi represents the charge fluctuations, 
while the current 3\ = ^dA\ describes the conserved fluc- 
tuations of vortices. As a consequence, the form B\ must 
be considered as a pseudo- vector, while A\ is a vector, as 
usual. The BF coupling is thus P- and T-invariant. 

The low-energy effective theory of the superconduc- 
tor can be entirely expressed in terms of the generalized 
gauge fields A\ and B\. The dominant term at long dis- 
tances is the BF term; the next terms in the derivative 
expansion of the effective theory are the kinetic terms 
for the two gauge fields (for simplicity of presentation we 
shall assume relativistic invariance), giving: 
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^dBi A *dBi , 



(4) 



where and are coupling constants of dimension m. 
Note that in (2-|-l) dimensions the action |4] is invariant 
under the duality transformation e g, A\ ■n- Bi. This 
action, including its non-Abelian generalization with ki- 
netic terms was first considered in [13 |. 



Naively one could diagonalize [4] by a transformation 
A = i(a -I- 5), B = (a-b), giving 

k f k f 

SBFid = 2) = — a A da bAdb. (5) 
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The result is a doubled Chern-Simons model for gauge 
fields of opposite chirality. It is the simplest example 
of the class of P- and T-invariant topological phases of 
strongly correlated (2+l)-dimensional electron systems 
considered in [T^ . 

The compactness of the gauge fields allows for the 
presence of topological defects, both electric and mag- 
netic. The electric (magnetic) topological defects couple 
to the form Ai (Bi) and are string-like objects described 
by a singular closed 1-form Qi (Mi). These forms rep- 
resent the singular parts of the field strenghts dAi and 
dBi, allowed by the compactness of the gauge symme- 
tries. The condensation of topological defects (or lack 
of) will lead to different topological phases characterized 
by a different ground state degeneracy on manifold with 
non-trivial topology. To show this we will use an ultra- 
violet (lattice) regularization. 

In we have shown that the zero-temperature par- 
tition function of (2-|-l)-dimensional BF may be written 
on the lattice as fields (vector) and (axial vector) 



Z = 



j^^^t^ j cxp(-5) 



{Qo) 

{Mo) 



dt^-l- A^K^^B^ 



-ikA^Qa + ikBoMo 



(6) 



Kfj^^ is the lattice Chern-Simons term [1^, defined by 
i^oo 0, Koi -£ijdj, Km = SiC^dj and K^j = 
—Sieijdo, in terms of forward (backward) shift and dif- 
ference operators Si (Si) and di (di). Its conjugate K^^^ 
is defined by Kqq = 0, Koi = —Sieijdj, Km = Cijdj and 
Kij = —SjCijOa. The two Chern-Simons kernels if^j/ 
and Kf^i, are interchanged upon integration (summation) 
by parts on the lattice. The topological excitations are 
described by the integer- valued fields Qo and Mq and rep- 
resent unit charges and vortices rendering the gauge field 
components Aq and Bq integers via the Poisson sum- 
mation formula; their fiuctuations determine the phase 
diagram The lattice spacing I is assumed ^ = 1. 

On the lattice the fields A^ and B^ are angular vari- 
ables [13] defined on the interval [— 7r,7r]. If we write 
A= i (o -I- 6) and i? = (a — 6) , it is now clear that the ul- 
traviolet regularization, required by the presence of topo- 
logical defects, make impossible to diagonalize naively 
the mixed Chern-Simons term. 

In the phase in which electric and magnetic topological 
defects condense the partition function requires a formal 
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sum also over the form Qi and Mi 
Z = VAVBVQVM 



k f 

expi— / (Ai A dBi + Ai A *Qi + Bi A *Mi)(7) 

Topological excitations may be absorbed into the com- 
pact gauge field A'^^, -B'^p. This case corresponds to the 
doubled Chern-Simons with both gauge fields compact 
described by the action [3] , with a ground state degener- 
acy (fci X k2)'^ on the torus. 

In case only one of the two topological defects con- 
denses we have two dual phases Q: when the magnetic 
excitations are dilute and the charge excitations condense 
rendering the system a superconductor: vortex confine- 
ment amounts here to the Meissner effect. If the magnetic 
excitations condense while the charged ones become di- 
lute: the system exhibits insulating behavior due to vor- 
tex superconductivity accompanied by a charge Meissner 
effect. Since the two phases are dual we will discuss only 
the charge condensation phase. 

First of all let us notice that associated with the con- 
finement of vortices there is a breakdown of the origi- 
nal U(l) matter symmetry under transformations Ai ^ 
Ai + dX. To see this let us consider the effect of such 
a transformation on the partition function ([T]) with an 
electric condensate. Upon integration by parts, the ex- 
ponential of the action acquires a multiplicative factor 



exp I 
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XA*Qi 



A A *Qj 



AIci,t— — oo 



(8) 

Assuming a constant A, we see that the only values for 
which the partition function remains invariant are 



A 
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n — 1 . . . ki 



(9) 



which shows that the global symmetry is broken from 
U(l) to Zk-^. In this phase all charges different from 
nkie are screened and the Dirac quantization condition 
imposes that diluted magnetic vortices are quantized in 
units r^n. Note that this is not the usual Landau mech- 
anism of spontaneous symmetry breaking. Indeed, there 
is no local order parameter and the order is characterized 
rather by the expectation value of non-local, topological 
operator the 't Hooft loop [HJ. The important point is 
that, in this phase, the ground-state degeneracy on the 
torus is not k'^ but rather only fc, as we will now show. 

By rewriting the charge topological excitations as the 
curl of an integer- valued axial field (p^: Qo = tijdiipj, the 
partition function become: 



TLB 



k f \ 

-i^ / dt22AoKoi{Bi + ^p^) 



From one sees that the gauge field components Bi 
are angular variables due to their invariance under time- 
independent integer shifts. Such shifts do not affect the 
last term in the action, which contains a time derivative, 
and may be reabsorbed in the topological excitations(pi, 
leaving also the first term of the action invariant. Topo- 
logical excitations may be absorbed only into the com- 
pact gauge B'^^ and one may define 

J2 J T^A^ j 'DB^ exp {^Sle {A„B,, ^,)) = 

J VA^,j VB% exp {^SLE{A^,B'i)) . (11) 

The canonical quantization of the low energy effec- 
tive action Sle {Ai^B'^i) is carried out by imposing the 
Gauss law constraints associated to the two Lagrange 
multipliers A^ and Bq, KoiB'^i = and KoiAi = and 
then enforcing the usual Weyl gauge condition = 0, 
B'^Q = 0: the canonical commutation relations read 
[Ai(x), Bj{x + i)] = — ^ij where i denotes a unit lat- 
tice vector in direction i. 

Since now only one of the two gauge fields is com- 
pact, we can introduce two pertinently normalized chi- 
ral gauge fields defined by A^i = {■^Ai + B'^i) and 
A^i = {^Ai — B'^i). Following [l^, the non-trivial wind- 
ings around closed loops can be encoded in the two pairs 
of global variables: 



k 

'2^ 



Y.A\{Q,l), 



^ — ' ZTT ^ — ' 



'A.KioBo + A,Ki.B. 



(10) 



where the sums run over a period Pi [l ~ 1,2) of the 
torus . As a consequence of the commutation relations 
between Ai and Bj one has that [qi,pj\ = iSij, i = R,L. 
The generators of large gauge transformations can thus 
be written as: 

Ui{n,m) = exp[2Tri{npL + mpi^)] 
U2{t,l) =exp[-ik{tqL + lqR)] , (13) 

and they satisfy the algebra: 

Ui{n, m)U2{t, I) — U2{t, l)Ui{n, m)xexp [i2TTk{nt + ml)] . 

(14) 

We note that for k an integer the cocycle is always trivial. 
Performing the transformation in the integer n,m,l, t: 

n = X + y ; m = x — y ] t = z + s \ I ~ z — s, (15) 

we can write 

Ui{n,m) = Ui{x)Ui{y) = 

= exp [2Trix{pL + Pr)] exp [2Triy{pL - Pr)] 

U2{t,l) = U2{s)U2{z)^ 

= exp [iks(qL + qR)] exp [ikz{qL - qR)] (16) 
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Ui{x) and Ui{y), 172 (s) and U2{z), are two sets of com- 
muting operators corresponding to the form and A^. 

Although one may define the winding numbers sepa- 
rately for the two chiral sectors, there is only one gen- 
erator of large gauge tansformations per homology cy- 
cle. This is because only the U{1) axial subgroup is com- 
pact; thus, large gauge transformations in the left sector, 
A^i A^i -\- diX with A {xi + P^) = X (xj) -I- 2nmi and 
rrii € Z for i — 1,2 must be combined with the corre- 
sponding large gauge transformations A^i — ^ A^^ + diX 
with X {xi + Pi) = X (xi) — 27rmi in the right sector. 

In [15] this implies x — z = making thus one of the 
generators in [16] act trivially on physical states. From 
here we can already see that the Hilbcrt space is only 
= fci X fc2 on the torus. We will however derive this 
explicitely for this gauge transformation. The generators 
of combined large gauge transformations (/l — 5- (7l + Stt, 
qR 9fl - 27r and pl ^ PL - k, PR Pn + k are 

C/i(l) =Ui = C/^iC/«i = exp [2m {pl - Pr)] , 
C/2(l) -C/2-C/^2C/'^2 = exp[zA:(9i-te)] , (17) 

and satisfy the algebra U1U2 = U2U1 exp [27rifc] . 

One is now ready to compute the ground state degen- 
eracy for a generic rational value of fc, i. c. k ~ ki/k2. 
On the torus, the wave function splits into a part depend- 
ing on the local variables (a;) and in a part depending on 
the global variables defined in ([T^ ; the degeneracy of the 
gro und state is entirely determined by this second part 
p^ . In the Schroedinger representation, the global vari- 
ables q are realized as generalized coordinates and the p 
are the corresponding momenta; thus, the wave function 
is a function of the q's only and must carry a representa- 
tion of the algebra of large gauge transformations (LGT). 
The representations of the algebra of LGTs are specified 
by two angles 61 and 62, which are entirely determined 
by the windings around the 2 homology cycles of the 
torus; correspondingly, the ground state wave function 
must also be labeled by these two angles, "0 = 1^6162- 
Assuming Ui ^pei62 — exp [i9i] ^'^162 one has 

U1U2 1^6162 = U2U1 exp [2Trikl] V'SiSa 

= exp [2nikl + iOi] Ui ^B,e2 , (18) 

and for I — k2 one obtains the same eigenvalue since 
exp [27rz(fci 7^2)^2] = 1, which leads to a first set of k2- 
degenerate independent states U2ipeie27 I = 1 • . . fe. For 
/ = fc2, one has U^'^ij^e^e^ = exp [102] i/'eiSa- 

The most general form of the wave function yielding a 
representation of LGT is 

i'ei92 = X! '^^P [*("^ + ^'i/27r)g] V'SiSa {n) , (19) 
nez 

where q — [q^ — qR)/2 is the normalized axial combina- 
tion of winding numbers. Using (|19p one has 

U2i>9ie2 = e:x.p[i{n + 0i/2Tr)q + ikql] V'SiSal'^) : 
nez 

(20) 



which reduces, for ^ = to 

exp [iq{n + 6'i/27r + ki)] ?/'eie2(«) = 

nez 

exp [26*2] '^exp[i{n + 0i/27T)q] ipe^e^ln) . (21) 
nez 

Combining this with [/2^Veie2 = exp [i6'2] ^'61^2 en- 
ables to derive the quasi-periodicity condition ijjeie2 (j^) — 
exp [1^2] '4'9i92i''^ + ki). For each I there are fci indepen- 
dent states; thus, the dimension of the Hilbert space is 
ki X k2 for the torus (generically {ki x k2)^ on genus g 
Riemann surfaces). 

For the superconducting and its dual insulator phase 
we have a degeneracy that is only ki x fc2, while when 
both gauge fields are compact the degeneracy is {ki x 
fe)^- From the phase structure analysis of BF topologi- 
cal fluids done in 0, 0| we see however that the phase in 
which both topological defects condense is not present, 
but, instead, we have either electric or magnetic con- 
densation or both are dilute. Thus the quantum order of 
the superconducting/insulator state BF topological fluids 
thus in the universality class deflned by the low-energy 
effective gauge theory 




-^Jd'x Ay^'^d^Al . (22) 

The action involves two separate Chern-Simons 

terms of opposite chirality. Contrary to previously 
considered examples of "doubled Chern-Simons theo- 
ries" [3], however, only the axial (vector) subgroup 
C^(l)axiai(vector) of the total gaugc group G = U{1)l ® 
U{1)r is compact; the diagonal vector group is the non- 
compact group R(L), reflecting the fact that the charge 
(vortex) coupled to it is no good quantum number due 
to the existence of a superconducting condensate. Thus, 
the corresponding topological order is halved: as a con- 
sequence BF topological fluids belongs to a universality 
class which is different than the one of conventional dou- 
bled Chern-Simons theories. 

We analyzed so far only the T = superconducting 
phase BF topological fluids; by means of duality one 
may derive the properties of the corresponding insulat- 
ing phase. What about the quantum transition point? 
To this end one should observe that a compact U{\) the- 
ory always involves a scale, determining the radius of the 
gauge group: thus, at the quantum critical point both 
gauge groups must decompactify and long-range quan- 
tum fluctuations are described by a continuum gauge 
theory with gauge group Rvector ® R-axiai- As a conse- 
quence, both charges and vortices must deconfine at the 
critical point, implying that, there, the BF topological 
fluids are- as expected - in a metallic phase; this be- 
havior supports the recent re-proposed scenario that 
quantum critical points are generically described by con- 
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tinuum gauge theories with deconfined degrees of free- dom. 
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